function y = A(t,a,b,x)
%   function y = A(t,a,b,x)
%
% t: valor (o vector de Mx1 de valores) en el rengo [0,1]
% a: vector en R^3x1
% b: vector en R^3x1
% x: matriz de tamanio 3xN que representa N valores en R^3 (en cada col)
% C: cota
% y: se retorna el valor de la funcion que indica el nivel de aturdimiento.

y=zeros(size(t),class(t));

for k=1:numel(t),
    
    ft = a * t(k) + b;
    resta = repmat(ft,1,size(x,2)) - x;
    norm2 = sqrt(sum(resta.*resta,1));
    y(k) = sum(1./norm2);

end

